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In low dimensional cuprates several interesting phenomena, including high T c superconductivity, are 
deeply connected to electron correlations on Cu and the presence of the Zhang-Rice (ZR) singlet state. Here, 
we report on direct spectroscopic observation of the ZR state responsible for the low-energy physical 
properties in two isostructural A-site ordered cuprate perovskites, CaCu 3 Co 4 0 12 and CaCu 3 Cr 4 0 12 as 
revealed by resonant soft x-ray absorption spectroscopy on the Cu L 3>2 - and O K-edges. These measurements 
reveal the signature of Cu in the high-energy 3+ (3d 8 ), the typical 2+ (3d 9 ), as well as features of the ZR 
singlet state (i.e., 3d 9 !, L denotes an oxygen hole). First principles GGA + U calculations affirm that the 
B-site cation controls the degree of Cu-O hybridization and, thus, the Cu valency. These findings introduce 
another avenue for the study and manipulation of cuprates, bypassing the complexities inherent to 
conventional chemical doping (i.e. disorder) that hinder the relevant physics. 

The many-body ZR singlet state 1 " 9 , a doped hole on the oxygen site coupled antiferromagnetically with a hole 
on the copper site, can be manipulated with conventional solid-state chemistry methods by partial removal 
of oxygen or by non-isovalent cation exchange (e.g. La 3+ —> Sr 2+ ). Such chemical routes of manipulation 
make the the doped cuprates prone to chemical disorder causing strong structural distortions or even changes in 
crystal symmetry, which in turn may severely alter the properties associated with Cu d-electron derived electronic 
and magnetic structures. En route to the goal of realizing the ZR state without chemical disorder or lattice 
distortion, while the majority of the Cu oxide compounds have the Cu ions in the B-site of the perovskite 
AB0 3 structure, the A-site ordered perovskites with chemical formula (ACu 3 )B 4 0 12 are intriguing candidates 
for investigation due to the unique Cu A-site arrangement (see Fig. 1(a)). In this structure, a surprisingly rich set 
of interesting physics phenomena 10 " 16 have been demonstrated by replacement of the B-site 3d transition metal 
ion via experiment and theory including colossal dielectric constant 1718 , ferrimagnetism 19 , charge, orbital, and 
spin ordering 20 , non-fermi liquid behavior 21 , enhanced electronic correlations 22 and possibly the Zhang-Rice 
state 23 . In addition, previous reports hinted at the presence of a mixed-valency on Cu for B = Cr and Co 23,24 . 
Despite the large variation in the B-site ions, all the members of this family belong to the same space group Im-3 
built upon two sub-lattices of octahedral B0 6 units and the planar Cu0 4 units (see Fig. l(b)-(c). As seen in Fig. l(a 
and d), the two structural units are connected via apical oxygens to form a B0 6 - Cu0 4 - B0 6 cluster. In this 
configuration, the effect of the B-site cation has been shown to be two-fold: first, its valence state controls the 
nominal valence state of the A' -site Cu ion and second, most importantly, its tendency for covalent mixing alters 
the charge distribution between Cu d- and O p-states 1115 . The structural similarity of the local ionic environment 
of Cu in these ordered perovskites with that of high T c cuprates, and, in close analogy to chemical doping, the 
above-mentioned B-site's potential to engineer the electronic structure prompts one to use the A-site ordered 
perovskites to gain unique insight into the development of the cuprate electronic structure. Construction of such 
manipulated valencies suggests the ability to tune the formation of ligand holes on oxygen and possibly to achieve 
the formation of the Zhang- Rice singlet state - an important ingredient for high temperature superconductivity. 
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Figure 1 | Crystal structure of A-site ordered perovskites. (a) The entire 
crystal structure of CaCu3Cr 4 0 12 . (b) Cr0 6 octahedral sub-lattice, 
(c) Cu0 4 distorted square planar sub-lattice, (d) Cr-Cr and Cr-Cu 
exchange pathways and (e) a sketch of B-site electronic states coupled to 
the ZR singlet state on Cu. 



Results 

Resonant x-ray absorption spectroscopy (XAS) in the soft X-ray 
regime is an ideal tool, extensively used in the past, to investigate 
the Cu and O electronic structure and the effects of electron -electron 
correlations in high T c cuprates 2,6 " 8,25 " 27 . Results of our soft XAS 
measurements on the Cu L 32 - and O K-edge are shown in Fig. 2 
for the two A-site ordered perovskites CaCu 3 Co 4 0 12 (CCCoO) and 
CaCu 3 Cr 4 Oi2 (CCCrO) along with that of the optimally doped 
superconducting YBa 2 Cu 3 0 7 - ( 5 (YBCO) and metallic LaCu0 3 
(LCO) perovskite with a unique formal Cu 3+ (3d 8 ) oxidation state. 
From the absorption line shape in Fig. 2, it is readily apparent that in 
CCCoO and CCCrO a mixed valency of Cu is present; we will discuss 
the CCCrO compound first. The sharp peak ~ 930 eV corresponds 
to the transition from the d 9 ground state, featuring a single hole in 
the Cu e^ band, to the cd 10 (here c denotes a core hole) excited state. 
The shoulder seen at 931.5 eV, a signature of the ZR singlet state, is a 
transition from the ground state d 9 L to the cd 10 L excited state 2,5 ' 7,26 ' 27 . 
The key difference between these two transitions being that while the 
Cu still maintains a nearly 2+ (or d 9 ) valency, a ligand hole is dis- 
tributed over the neighboring oxygens as illustrated in Fig. 1(e). The 
interaction between this ligand hole and the core hole created by the 
photon absorption raises the energy required to promote the core 
electron to the unoccupied state in the e^ band resulting in the 
observed high- energy state marked by red arrows. Furthermore, to 
corroborate this, we show the absorption data on LCO with the 
formal Cu charge state of 3 + . Indeed, in-line with our previous 
statement, the main feature of the LCO XAS around the L 3 -edge at 
— 931.5 eV is the transition from the d 9 L state, thus making the ZR 
state the dominant contribution to the ground state in good agree- 
ment with previous reports 28 ' 29 . The aforementioned doped ligand 
holes are the charge carriers essential to the physics of high T c cup- 
rates; for instance, in hole doped cuprates the superconducting trans- 
ition temperature is strongly dependent on the amount of doping 
present 2,30,31 . To clarify the connection between these compounds 
and the high-temperature superconducting cuprates, XAS data was 
taken for an optimally doped YBCO sample. As seen in Fig. 2, the Cu 
L-edge spectrum of CCCrO shows remarkable resemblance to 
YBCO, namely the dominant d 9 initial states followed by the higher 
energy ZR state shoulder. This result highlights the capacity of the A- 
site ordered cup rate perovskites to mimic the effects of chemical dop- 
ing without provoking unwanted lattice and stoichiometry devia- 
tions inherent to conventional chemical doping. In addition, the 
multiplet split peak at ~ 940 eV corresponds to a transition from 
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Figure 2 | XAS measurements on Cu L- and O K-edges. Left panel 
displays the Cu L-edge for all samples with the Zhang-Rice and d 8 ground 
states highlighted. LCO and YBCO are shown as references for Cu 3 + and 
Cu 2+ respectively, highlighting the change in valence for Cu with B-site 
change 36,37 . Right panel shows the corresponding O K-edge spectra 
displaying the Zhang-Rice ligand hole state. In LCO, the higher energy 
peak around ~ 940 eV is due to the presence of the unusual Cu d 8 state. In 
CCCoO, a previously reported impurity peak around 929 eV 32 is 
subtracted from both the L 3 and L 2 edges by fitting to a Voigt function. This 
peak corresponds to the d 9 peak observed on the other samples, however it 
has a very small spectral weight (~ 1/50 ratio between d 9 and d 9 L states). 

the metastable 3d 8 state to the cd 9 excited state, indicating the pres- 
ence of the ionically Cu 3+ state entirely absent in YBCO. 

Next, we turn our attention to the CCCoO electronic structure. As 
clearly seen in Fig. 2, CCCoO shows a dominant white absorption 
line analogous to CCCrO, but, surprisingly, it is shifted ~ 1 eV away 
from the d 9 white line of CCCrO and YBCO. When compared with 
the formally Cu 3+ compound LCO (shown immediately above it in 
Fig. 2) it becomes readily apparent that the Cu ground state is pre- 
dominantly of d 9 L character, implying a dominant Cu 3+ formal val- 
ence as in LCO. To quantify the prevalence of the d 9 L state, a three 
peak fit is used for the L 3 data, which are centered at 929.3 eV 
(impurity) 32 , 930 eV (d 9 ), and 931.4 eV (d 9 L) which gave a very small 

spectral weight ratio of .0203 (^j^j • The impurity peak was sub- 
tracted from the original data (dotted line) for clarity. Further dem- 
onstrating the 3+ formal valence, the multiplet split peak at 
—940 eV provides strong evidence for the admixture of the d 8 state. 
When considered alongside CCCrO, CCCoO also contains all three 
Cu states, but with obviously significantly different weights, further 
testifying to the ability of the B-site to effectively hole dope the Cu 
sites. 

Because of the large degree of charge transfer between oxygen and 
copper and its importance to the ZR state, studies of the O K-edge 
can lend supporting information about the electronic structure of 
Cu. As seen in Fig. 2(b), O K-edge data in all samples show the 
oxygen specific signature of the 3d 9 L state, a necessary ingredient 
of the ZR state, which is particularly evident from the comparison to 
the well studied case of YBCO, where the ligand hole on oxygen 
appears as a large, doping dependent pre-peak around 528 eV 6 " 8,25 . 
This observation of the pre-peak intensity indicates that the ZR state 
is the first ionization state at the Fermi energy. This further evidences 
our conclusion on the presence of the ZR state inferred from the Cu 
L-edge analysis. 

Although our prime focus is on the electronic structure of Cu and 
O, Co and Cr L 3j2 -edge absorption can serve to independently 
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Figure 3 | XAS measurements on the Cr and Co L-edges. (a) The Cr L- 

edge along with 3+ and 4+ references from 33 , the CCCrO and the 4+ 
reference samples show strong similarities, (b) The Co L-edge with valence 
references from 34 showing a similar, but not identical line shape. No 
standards were available meaning any energy shift between the samples 
could not be analyzed, all spectra were shifted to ~ align the white lines. 

validate the Cu oxidation state. With this goal in mind, we performed 
L-edge XAS measurements on Co and Cr and compared this to the 
previous results on reference Cr and Co compounds. As shown in 
Fig. 3(a), a direct comparison to the Cr 3+ and Cr 4+ absorption con- 
firms that in CCCrO the Cr valency is 4+ or very close to it 33 . This 
verifies the conclusion that, in this compound, if Cr is 4+ then Cu 
must be —2 + (assuming O 2- and Ca 2+ ). Any deviation of Cr valency 
towards 3+ must then be compensated by raising the Cu charge 
state, thus Cu must be 2 + or very close to it, consistent with the 
Cu L-edge results. The Co L-edge data shown in Fig. 3(b) bear strong 
similarities with previous studies of LaCo0 3 , which consists of ion- 
ically bound Co 3+34 . Specifically, strong, high energy multiplet 
shoulders on the L 3 and L 2 white lines are present in both CCCoO 
and LaCo0 3 . In variance to this, Co 4+ shows a smaller shoulder 
shifted to even higher energy. While the high energy side of the 
CCCoO Co L-edge is consistent with Co 3+ , the low energy multiplet 
peak does not match either 3+ or 4+, implying CCCoO is neither 
ionically pure 3 + or 4+ . Based upon the Cu L-edge conclusion of the 
Cu 3+ charge state, we can independently assess the Co valency to be 
~ 3.25 + . In this context, the Co L-edge data is compatible with the 
Cu L-edge results and implies the presence of a mixed 3 + /4+ Co 
valency, although it cannot confirm the exact 3.25+ valence state 
needed to conserve charge. Lending further support to this conclu- 
sion, Co L edge data in total electron yield (TEY) mode is shown in 
the Supplemental Figure SI. 

To gain further understanding of the B-site influence on the 
unusual electronic states on Cu, we performed first-principles 
spin-polarized DFT calculations within the GGA+ U approximation. 
Fig. 4 summarizes the first-principles results for CCCoO and 
CCCrO, respectively. Note, although no long-range magnetic order- 
ing is reported in these compounds, spin-polarized calculations are 
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Figure 4 | The GGA + U density of states, projected onto Co d (Cr d), Cu 
d and O p states for the A-site ordered perovskites. (a) CaCu 3 Co 4 0 12 
((b) CaCu3Cr 4 0 12 ) compound. The positive and negative values in y-axis 
correspond to density of states inup and down spin channels, respectively. 
The zero of the energy is set at E F . The features due to dominant 
contributions of various degrees of freedom have been marked. The inset 
in the middle panel shows the plot of the effective Cu x 2 — y 2 Wannier 
function, calculated by NMTO-downfolding. Shown are the isosurface 
with lobes of opposite signs are colored differently. The inset in the 
bottom panel shows the relative positions of Co d, Cu d and O p energy 
levels. 

essential for proper assignment of both charge and spin states of 
various constituents 15 . The three panel Fig. 4(a) and (b) show an 
inspection of the calculated density of states (DOS) projected onto B- 
d (B = Co/Cr), Cu-d and O-p states respectively. A direct inspection 
of the states in the vicinity to the Fermi energy, E F demonstrates the 
existence of strong mixing between B and Cu d-states and O p-states. 
This strong covalency effect makes both CCCoO and CCCrO sys- 
tems metallic, in agreement with experimental observations. In addi- 
tion, as shown in the inset of Fig. 4(a) -(b), the octahedral crystal field 



SCIENTIFIC REPORTS | 3 : 1834 | DOI: 1 0.1 038/srep01 834 



3 



of the oxygen atoms around the B ion splits the B-t 2g states from the 

B-e^ (e°^J states; the trigonal distortion present in the B0 6 octahedra 

breaks the t 2g levels further into doubly degenerate, e^ and singly 
degenerate a. lg states. The distorted square planar environment of 
oxygen atoms surrounding the Cu ion, on the other hand, lifts the 5- 
fold degeneracy of Cu d-levels, with large energy separation of Cu 
d x 2_ y i state from the rest. Here we note that the nominal d 8 (i.e. 
Cu 3+ ) would imply empty Cu: d x 2_ y 2 states in both spin up and spin 
down channels, and for d 9 Cu : d x 2 _ y i states would be occupied in 
one of the spin channels and be empty in the other spin channel. 
From the DOS plots shown in Fig. 4(a) and (b), we find that Cu: 
d x 2 _ y i states are empty in the down spin channel, while in the up spin 
channel it is indeed close to being empty for CCCoO, but mostly 
occupied for CCCrO. We note that DFT calculation deals with single 
configuration nature of the wavefunction and therefore, delivers the 
state with most dominant weight. Interestingly enough, the dom- 
inance of Cu 3+ state in CCCoO and Cu 2+ state in CCCrO is 
evident even at the level of single configuration theory. This result 
lends strong theoretical support for the conclusion inferred from our 
XAS measurements that the Cu state is predominantly in the 3 + 
state for CCCoO and closer to 2+ for CCCrO. 

Next, we discuss the calculated magnetic states in connection to 
Cu. First, we note that occupancy of the square planar Cu d-states 
according to the formal valence of d 8 would give rise to a paired 
electronic configuration with zero magnetic moment at Cu sites, 
while that of formal valence of d 9 would give rise to one unpaired 
electron with magnetic moment of 1 \i B in the fully spin polarized 
limit. In comparison, the calculation yields magnetic moments on Cu 
site 0.07 }i B for CCCoO and 0.50 /n B for CCCrO, thus lending addi- 
tional validation of our experimental conclusion on the exotic 3 + 
state for CCCoO and closer to 2+ valency in CCCrO. Moreover, 
from the simple ionic charge count, the nominal 3 + and 2 + valences 
of Cu (assuming Ca 2+ and O 2 "), would set the formal valences of the 
B-site ion to 3.25+ and 4+ respectively, and result in d-electron 
occupancies of 5.75 for Co and 2 for Cr. The computed magnetic 
moment at the Cr site is 2.2 /ig, in good agreement with the nearly 
2+ valence of Cu and 4+ valence of Cr in CCCrO. The d- occupancy 
of 5.75 on Co can lead to either low or high spin magnetic state; our 
computed Co magnetic moment of —1.7 /ig, however, is in confor- 
mity with the intermediate rather than low- or high-spin state. 

By examining the energy level positions, as calculated in DFT, 
shown in the insets in the bottom panels of Figs. 4(a-b), we find 
that Cu d x 2_ y i energy level is nearly degenerate with O-p states, 
leading to a situation very similar to that in high T c cup rates. This 
electronic configuration naturally supports strong mixing of Cu 
d x i_ y 2 and O p x and p y states by forming a strong pda anti- 
bonding combination, which may further bind to d-states at the 
B site as shown in Fig. 1(e). The formation of such a mixed state is 
also evident in the Wannier function plot of an effective Cu x 2 — 
y 2 Wannier function obtained by keeping only Cu d x 2_ y 2 state 
active and downfolding or integrating out all the other degrees 
of freedom (see inset of the middle panel of Fig. 4). We notice that 
for the Co compound mixing of Cu-O hybridized states with B- 
site d-states is much stronger compared to that of Cr compound, 
which is driven by the strong hybridization of empty C\x x 2_ y 2 
orbital of Cu 3+ state with empty Co, e^ orbitals. In CCCrO, this 
leads to vanishing contribution of the B-site states for the selected 
iso-surface value of the Wannier plot, showing the typical Cu-O 
antibonding wavefunction, making the situation more like high T c 
cuprates with dominant Cu 2+ state. The strong mixing between 
C\x x 2_ y 2 and O p states, which is a prerequisite for formation of 
the ZR state - observed as the feature around 931 eV in Cu L-edge 
spectra, is further corroborated by the presence of a 0.2 fi B frac- 
tion of the magnetic moment on the O site and 0.5 /ig on Cu. We 
thus conclude that the formation of a ZR like many-body state 



similar to what is observed on the high T c cuprates is indeed 
highly favorable in these systems with enhanced probability for 
the Cr compound. 

Discussion 

In summary, resonant soft X-ray absorption data on Cu L 3)2 - and O 
K- edges revealed the presence of the ZR singlet state in the symmet- 
rically invariant A- site- ordered cuprate compounds, CaCu 3 Co 4 0 12 
and CaCu 3 Cr 4 0 12 . Additionally, the data demonstrated the existence 
of a mixed valency for Cu modulated by the electronic state of the B- 
site cation and reflected in the d 9 , d 9 L, and the metastable d 8 states of 
Cu d- electrons. First principles GGA + U calculations further val- 
idate the notion that within the A-site ordered perovskite family, the 
choice of B-site cation is an effective tool to manipulate the degree of 
B-Cu-O hybridization, adjusting the Cu valency towards the unusual 
3 + valence state in CCCoO and away from 3 + closer to the stable 
2+ valency in CCCrO. Additionally, Wannier function calculations 
support this direct experimental observation that in the CCCrO 
compound the ZR many-body state is particularly favorable and 
controls the low-energy physics. These findings demonstrate that 
the charge and spin state of Cu, fundamental to the intriguing phys- 
ical properties the cuprates display, can be effectively altered by a 
careful choice of the B-site cation on isomorphic lattices that circum- 
vent the chemical disorder and lattice modulation intrinsic to other 
doping methods. 

Methods 

Experimental details. Polycrystalline CaCu 3 B 4 Oi2 (B = Cr and Co) samples were 
prepared under high-pressure and high-temperature (HPHT) conditions with a 
Walker-type multianvil module (Rockland Research Co.). For B = Cr, stoichiometric 
amounts of CaO, CuO, and Cr0 2 powders were thoroughly mixed and then subjected 
to a HPHT treatment at 7 GPa and 900°C for 30 minutes. For B = Co, a precursor 
containing Ca, Cu, and Co ions in the 1:3:4 molar ratio obtained via a sol-gel route 
was mixed with 30 wt.% KC10 4 acting as an oxidizing agent and then subjected to a 
HPHT treatment at 9 GPa and 1000°C for 30 minutes. The resultant KC1 was washed 
away with deionized water. Details about the sample assembly and procedures for the 
HPHT experiment can be found in a previous paper 35 . Phase purity of the above 
samples was examined with powder X-ray diffraction (XRD) at room temperature 
with a Philips Xpert diffractometer (Cu Ka radiation). The B = Co sample contains a 
small amount of cobalt and copper oxides. The cubic lattice parameter was calculated 
to be a = 7.2385(5) and 7.1259(3) A for M = Cr and Co, respectively. Reitveld 
refinement for the two a-site ordered compounds is available in the Supplemental 
Figure S2. 

LCO was synthesized under high oxygen pressure to the rhombohedral symmetry, 
details of which can be found elsewhere 36 . YBCO was grown with pulsed laser 
deposition, the details of which can also be found elsewhere 37 . 

In order to elucidate the transformation of the Cu valency across these compounds, 
soft XAS measurements were carried out at the soft x-ray branch at the 4-ID-C 
beamline in the bulk- sensitive total fluorescence yield (TFY) mode and TEY mode in 
the Advanced Photon Source at Argonne National Laboratory. Measurements were 
taken on the Cu L-edge and O-K edge for all samples. All measurements were 
accompanied by a CuO standard, which allows precise alignment of the energy for 
oxidation state comparison. The Co and Cr L-edges were also obtained for the 
corresponding samples. All measurements shown were obtained in TFY mode; 
except, the Cr L edge data and all the Cr reference absorption data were taken in TEY 
mode. 

Computational details. In the first-principles DFT calculations we have used the 
plane wave basis set and pseudo-potentials as implemented in the Vienna Ab-initio 
Simulation Package (VASP) 38 . The exchange-correlation function was chosen to be 
that of generalized gradient approximation (GGA) implemented following the 
parametrization of Perdew-Burke-ErnzerhoP 9 . The electron-electron correlation 
beyond GGA was taken into account through improved treatment of GGA + U 
calculation within the + U implementation of Dudarev et al. 40 . For the plane wave 
based calculations, we used projector augmented wave (PAW) 41 potentials. The wave 
functions were expanded in the plane wave basis with a kinetic energy cutoff of 
600 eV and Brillouin zone summations were carried out with a 6 X 6 X 6 k-mesh. 
The U value of 5 eV on Cu site and 4 eV on B site were used for GGA + U calculations 
while the Hund's rule coupling / was fixed to 0.8 eV. The obtained results were 
verified in terms of variation of U parameter. 

In order to estimate the positions of the Cu-d, B-d and O-p energy levels as 
well as the plots of the effective Wannier functions for B-d states, we used muffin-tin 
orbital (MTO) based N-th order MTO (NMTO) 42 -downfolding calculations. Starting 
from a full DFT calculations, NMTO- downfolding arrives at a few- orbital 
Hamiltonian by integrating out degrees which are not of interest. It does so by 
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defining energy- selected, effective orbitals which serve as Wannier-like orbitals 
defining the few-orbital Hamiltonian in the downfolded representation. NMTO 
technique which is not yet available in its self-consistent form relies on the 
self- consistent potential parameters obtained out of linear muffin-tine orbital 
(LMTO) 43 calculations. The results were cross-checked among the plane wave and 
LMTO calculations in terms of total energy differences, density of states and band 
structures. 
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